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INTRODUCTION 


During  the  past  few  years,  it  has  become  apparent  that 
problems  of  wave  propagation  in  the  earth  cannot  be  adequately 
studied  without  appropriate  information  regarding  the 
mechanical  behavior  of  the  in  situ  material.  Refs.  [1],  [2], 

and  [3].  Any  or  all  of  the  practices  of  co r ing,  samp  1 i n g , 
sample  preparation  and  testing  may  introduce  substantial  un- 
certainties and  errors  into  any  laboratory  assessment  of  the 
behavior  of  a particular  geological  site.  For  this  reason, 
various  in  situ  testing  procedures  have  been  developed  to 
directly  determine  the  ground  response  to  various  types  of  loading. 

This  data  can  be  used  in  conjunction  with  laboratory  data  to 
obtain  a more  complete  picture  of  the  in  situ  material  behavior. 

One  of  the  more  promising  of  the  dynamic  in  situ  tests  is 
the  Cylindrical  In  Situ  Test  (CIST),  Ref.  [4].  In  such  a test, 
a vertical  cylindrical  cavity  two  feet  in  diameter  is  explosively 
loaded  (and  subsequently  unloaded  as  the  cavity  pressure  decays) 
in  an  axisymmetric  fashion.  The  cavity  generally  extends 
through  several  layers  of  geological  materials  to  depths  ranging 
from  30  to  75  feet.  A number  of  vertical  boreholes  at  various 
ranges  and  azimuths  from  the  central  cavity  contain  instrumenta- 
tion at  a number  of  depths  to  measure  ground  motion  in  various 
directions  in  the  different  layers  of  material. 

The  purpose  of  the  present  study  is  to  investigate  the 
means  by  which  the  results  obtained  in  CIST  experiments  can 
be  utilized  in  the  construction  of  a material  model  to  rep- 
resent in  situ  behavior.  Because  the  CIST  configuration  does 
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not  exercise  the  material  in  all  stress  and/or  strain  paths 
expected  in  a ground  shock  problem,  the  use  of  data  from 
laboratory  tests  is  also  incorporated  into  the  modeling  pro- 
cedure as  described  in  this  report. 

The  cap  model.  Refs.  [ 5 ] - [ 7 ] , was  used  for  this  study 
because  many  recent  ground  shock  computations  have  utilized 
this  type  of  model.  The  model  is  capable  of  representing  a 
wide  range  of  material  behavior  while  satisfying  conditions 
of  existence,  uniqueness  and  stability  of  solution.  In  the 
present  study,  the  capability  of  the  cap  model  to  represent 
material  behavior  in  CIST  is  investigated. 


1 1 OUTLINE  OF  STUDY 

The  procedure  employed  in  this  study  is  based  on  the 
premise  that,  if  possible,  the  material  model  should  represent 
as  many  important  features  of  material  behavior  as  is  practi- 
cal. Specifically,  one  should  select  the  manner  in  which  the 
CIST  test  is  employed  for  the  construction  of  a model  in  such 
a way  as  to  minimize  the  differences  between  the  behavior  of 
the  resulting  model  and  the  laboratory  based  property  data 
(provided,  of  course,  there  is  no  inherent  inconsistency  be- 
tween the  CIST  and  laboratory  results). 

The  most  natural  way  to  accomplish  this  objective  is  to 
define  the  qualitative  behavior  of  the  material  in  a number 
loading  and  unloading  paths.  Because  materials  can  be  subjected 
to  a wide  range  of  controlled  stress  and/or  strain  paths  in 
the  laboratory,  the  laboratory  data  will  lead  to  a selection 
of  model  functional  forms  and  parameters  which  properly  represent 
the  variations  in  material  behavior  under  these  paths.  The 
CIST  measurements  can  then  be  used  to  adjust  the  model  parameters 
(if  possible)  to  arrive  at  a final  model  which  is  in  quanti- 
tative agreement  with  in  situ  behavior. 

The  starting  point  for  the  present  study  involves  fitting 
the  cap  model  to  available  laboratory  data  for  the  site  under 
consideration  - the  HARD  PAN  site  in  Kansas  at  which  CIST  12 
was  performed.  Using  the  laboratory  based  cap  model,  a cal- 
culation of  the  CIST  event  was  performed  and  the  results 
compared  to  the  measurements  made  during  CIST.  Based  on  these 
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comparisons,  several  changes  were  made  In  the  model  parameters 
and  the  calculations  were  rerun.  An  Iterative  procedure 
involving  computation,  comparison  of  results  and  measurements, 
and  alteration  of  the  model,  was  undertaken  in  order  to 
reduce  the  discrepancies  between  experiment  and  calculation 
until  satisfactory  CIST  results  were  obtained  in  the  computation. 
As  expected  it  was  necessary  to  relinquish  some  of  the  quanti- 
tative agreement  of  the  model  and  laboratory  data  in  order  to 
achieve  improved  results  in  the  CIST  computation.  The 
differences  between  the  final  CIST  based  model  and  the  labora- 
tory data  are  not  unduly  large,  however,  and  are  considered 
to  be  well  within  the  usual  uncertainties  involved  in  representing 
material  behavior  at  a real  site.  The  CIST  based  cap  model 
obtained  in  this  study  is  currently  being  used  to  represent 
the  HARD  PAN  site  in  free-field  and  structure-medium  interaction 
problems  for  the  Space  and  Missile  Systems  Organization  (SAMSO). 


Ill  SITE  PROFILE  AND  MODEL 


The  particular  site  chosen  for  this  study  is  the  HARD  PAN 
site  in  Kansas  at  which  CIST  12  was  performed.  The  HARD  PAN 
series  of  field  tests  is  designed  to  investigate  the  response 
of  a scaled  Minuteman  Launch  Facility  to  explosively  generated 
air  blast  and  ground  shock  in  an  inter  bedded  clay/shale/limestone 
geology.  The  HARD  PAN  geology  is  complex  and  varies  significantly 
from  point  to  point  over  the  site.  An  idealized  and  simplified 
site  profile  based  on  Refs.  [8]- [10]  was  chosen  for  the  CIST 
calculations  performed  in  this  study  in  an  attempt  to  "average" 
the  material  properties  for  each  layer  in  which  the  CIST 
measurements  were  taken.  There  are  two  reasons  to  make  such  a 
simplification.  First,  the  exact  site  profile  is  not  known  for 
each  CIST  guage  hole;  it  varies  somewhat  with  aziumth  and 
range  from  the  central  CIST  cavity.  Second,  two  dimensional 
axisymmetric  calculations  of  practical  size  cannot  adequately 
represent  that  portion  of  the  ground  signal  arising  from  such 
small  scale  variations  in  material  properties.  The  idealized 
HARD  PAN  site  profile  chosen  for  this  study  is  shown  in  Fig.  1. 

A set  of  laboratory  based  data  which  describes  the  material 
behavior  of  each  of  the  layers  of  Fig.  1 was  taken  from  Ref.  [8]. 
Although  a more  detailed  set  of  material  properties  is  included 
in  Ref.  [9],  the  use  of  such  a set  of  properties  for  the  present 
study  is  not  warranted  for  the  reasons  described  above.  The 
laboratory  based  material  behavior  for  each  of  the  five  materials 
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LAYER  1 

DRY  CLAY  < 

LAYER  2 

WET  CLAY 

LAYER  3 

LIMESTONE 

LAYER  4 

SHALE 

LAYER  S 

COAL 

LAYER  6 

SHALE 

LAVER  7 

LIMESTONE 

idealized  site  profile  for  hard  pan 


included  in  the  current  study  was  fit  by  means  of  an  eight 
parameter  version  of  the  cap  model.  Ref.  [5],  modified,  where 
required,  to  represent  rocks  as  described  in  Refs.  [6]  and  ( 7 ] . 


The  model,  which  is  depicted  in  Fig.  2,  includes  a failure 
envelope, 

Sj'2  = FF(J1)  = A - C exp(BJj)  (1) 

an  elliptical  cap  of  constant  eccentricity, 

/j'2  = Fc(J1,k)  = | Z[X(k)-L(k)  ]2-  [ J j^-L  (k)  J2  (2) 

which  is  governed  by  the  hardening  rule 


ej  = W [ exp  (DX)  - 1]  (3) 

and  constant  bulk  and  shear  moduli,  K and  G , within  the  yield 
surface.  In  the  above  equations  A,  B,  C,  R,  W and  T)  are  model 
parameters,  J ^ is  the  first  invariant  of  stress  (tension 

f 

positive),  J2  is  the  second  invariant  of  the  stress  deviators, 

K is  the  cap  hardening  parameter  and  X and  L are  functions  of 
K through  the  relations 


L 


[ 


K if  K < 0 

0 if  k > 0 


X - K - R[A  - C exp(B<) ) 


(A) 


The  quantity  depends  on  the  plastic  volumetric  strain,  , 

through  the  relations 
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FIG.  2 YIELD  SURFACES  IN  CAP  MODEL 


In  this  study  the  clays,  i . e . , the  two  upper  layers  in  Fig.l, 
are  considered  to  be  soils.  The  material  is  assumed  to  be 
incapable  of  supporting  hydrostatic  tension  - all  stresses 
components  are  set  to  zero  when  hydrostatic  tension  occurs. 
The  model  is  described  in  greater  detail  in  Ref.  [7]. 

The  parameter  values  for  which  the  cap  model  fits  the 
laboratory  based  material  behavior  are  listed  in  Table  I,  and 
the  behavior  of  this  model  is  shown  in  Figs.  (3-12)  together 
with  the  laboratory  data  to  which  the  model  was  fitted. 


TABLE  1 

CAP  MODEL  BASED  ON  LABORATORY  DATA 

MATERIAL 


MODEL 

DRY 

WET 

" ~ 

PARAMETER 

CLAY 

CLAY 

LIMESTONE 

SHALE 

COAL 

K (ksi  ) 

723. 

725  . 

1 600  . 

1160. 

290. 

' G(ksi) 

2 . 9 

4.35 

957  . 

24  7 . 

174. 

A ( k s 1 ) 

. 0334 

.016 

58. 

55.2 

33 . 8 

B(ksi_1) 

0. 

0. 

. 00966 

. 00966 

. 00966 

C ( k s 1 ) 

0. 

0. 

57.  1 

55.  1 

33.4 

D(ksl_1 ) 

.483 

3.43 

.69 

.11 

.586 

W 

.054 

.0021 

.00025 

.0029 

. 002 

R 

2.3 

2 . 5 

3.5 

5.0 

3 . 8 

Density  (p.c.f.) 

124. 

128. 

167  . 

155  . 

84 . 9 

I 
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OTX)  (ISX) 


UNIAXIAL  STRAIN  TEST  HARD  PAN  LIMESTONE 

►—  LA1-MDDEL 
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FIG. 7 COMPARISON  OF  LAB  DATA  AND  LAB  BASED  MODEL 
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UNIAX  STRESS  PATH  HARD  PAN  LIMESTONE 

* LA1-NBDCL 

LAB • DATA 
FAILURE 


PRESSURE  (KSI) 


FIG.  8 COMPARISON  OF  LAB  DATA  AND  LAB  BASED  MODEL 
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UNIAXIAL  STRAIN  TEST  HARD  PAN  SHALE 


* LAS-MKL 

lm<  data 
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FIG.  9 COMPARISON  OF  LAB  DATA  AND  LAB  BASED  MODEL 
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FIG.  10  COMPARISON  OF  LAB  DATA  AND  LAB  BASED  MODEL 
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AXIAL  STRESS  (KSI) 


UNIAXIAL  STRAIN  TEST  HARD  PAN  C0AL 


COMPARISON  OF  LAB  DATA  AND  LAB  BASED  MODEL 
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FIG.  II 
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FIG.  12  COMPARISON  OF  LAB  DATA  AND  LAB  BASED  MODEL 
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IV  COMPUTATIONAL  SETUP 

The  model  described  in  the  proceeding  section  was  used 
to  perform  a calculation  of  the  CIST  Event.  For  this  purpose 
the  LAYER  code,  which  is  briefly  described  in  the  Appendix, 
was  used.  The  computational  grid  layout  is  shown  in  Fig.  (11). 
The  time  step  for  the  calculation  was  0.08  ms  and  the  calcula- 
tion was  carried  to  a time  of  about  10  ms.  The  computation 
required  approximately  six  minutes  of  central  processor  time 
on  a CDC  8600. 

The  proper  input  loading  for  the  CIST  computation  is  a 
matter  of  considerable  uncertainty.  No  reliable  pressure 
measurement  was  obtained  for  the  HARD  PAN  CIST  test  and 
considerable  variation  in  the  cavity  pressure  loadings  has  been 
observed  in  other  CIST  tests.  Therefore  it  was  decided  that 
the  nominal  or  idealized  cavity  pressure  loading  shown  in 
Fig.  14  would  be  used  unless  it  was  later  found  to  be  incorrect 
or  inadequate.  During  the  present  study  no  alteration  of  the 
input  loading  of  Fig.  14  was  found  to  be  required. 

Several  checks  were  performed  to  validate  the  LAYER  code 
setup  to  handle  CIST.  One  and  two  dimensional  axisyrametrlc 
calculations  were  performed  with  the  bilinear  material  model 
of  Ref.  [11].  In  these  calculations  variations  were  made  in 
the  grid  size  and  the  procedure  used  to  apply  the  cavity 
pressure  boundary  condition  to  the  grid.  The  results  of  these 
computations  were  compared  to  results  obtained  by  means  of  the 
characteristic-based  code  of  Ref.  [11).  The  agreement  between 
the  characteristic  code  and  LAYER  was  quite  good,  as  reported 
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in  Ref.  [111. 


An  additional  check  on  the  effect  of  grid  size  was  made 
by  rerunning  the  early  portion  of  the  cap  model  CIST  calcula- 
tion with  a 1/2  foot  radial  mesh  spacing  instead  of  the  1 foot 
spacing  originally  used.  No  significant  differences  were 
observed  in  the  results,  so  the  original  1 foot  square  mesh 
was  retained  for  the  entire  study  performed  here. 

As  indicated  in  the  Appendix,  LAYER  is  a small  displacement 
code,  i.e.,  the  assumption  is  made  in  this  study  that  convection 
and  transport  effects  near  the  CIST  cavity  can  be  neglected. 
Furthermore,  the  small  air  blast  loading  effects  on  the  ground 
surface  away  from  the  cavity  are  also  neglected,  and  the  applied 
pressure.  Fig.  14,  is  assumed  to  fce  uniform  throughout  the 
cavity.  The  agreement  obtained  between  the  computations  and 
the  measured  CIST  data  - to  be  shown  in  the  next  section  - 
indicates  that  the  above  assumptions  are  appropriate  as  well 
as  extremely  convenient.  This  is  important  because  these 
assumptions  drastically  reduce  the  computer  time  required  to 
perform  a CIST  computation  and  therefore  enable  one  to  perform 
the  iteration  procedure  required  in  this  study  quickly  and 
efficiently. 
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V 


THK  ITERATION  PROCEDURE  AND  CIST  BASED  MODEL 


The  iteration  procedure  used  in  this  study  to  develop 
the  final  ('.1ST  based  model  consists  of  the  following  three 
steps  : 

1.  Compare  the  results  of  I he  CIST  calculations  with 
measured  CIST  ground  response,  noting  similarities 
and  differences  in  the  waveforms  as  well  as  the 
sensitivity  of  various  waveform  features  to  model 
changes  made  in  previous  CIST  calculations. 

2 . Based  on  the  above  comparisons,  determine  physically 

acceptable  model  parameter  changes  which  are  likely  to 
reduce  the  discrepancies  between  the  CIST  calculations 
and  the  CIST  experimental  data. 

3.  Make  the  required  revisions  in  the  material  model 
and  recalculate  the  CIST  event. 

These  three  steps  are  repeated  until  either  a)  satisfactory 
agreement  between  calculation  and  experiment  is  obtained,  or 
b)  Inability  of  the  model  to  represent  the  material  behavior 
in  the  CIST  event  is  demonstrated. 

The  iterative  procedure  employed  in  this  study  begins  with 
the  comparison  of  the  CIST  measurements  and  the  results  of  the 
CIST  calculation  performed  w i L li  the  lab  based  model.  A typical 
set  of  waveforms  used  for  this  comparison  is  shown  in  Fig.  (15), 
in  which  the  horizontal  particle  velocity  versus  time  is  plotted 
at  ranges  of  3 feet,  5 feet  and  8 feet  at  a depth  of  45  ft. 


(in  the  coal  layer).  The  computed  waveforms  are  qualitatively 
quite  similar  to  the  measured  CIST  waveforms,  especially  if 
one  allows  for  baseline  shift  corrections  to  the  CIST  data. 

The  major  discrepancies  between  the  calculated  and  measured 
waveforms  occur  in  the  magnitude  of  the  peak  velocity,  which 
is  too  high  in  the  computations,  and  the  arrival  time  of  the 
wave  front,  which  is  too  late  in  the  computations. 

These  facts  suggest  that  the  behavior  of  the  lab  based 
model  is  too  soft  for  initial  loading  (a  stiffer  loading  modulus 
tends  to  reduce  peak  particle  velocity  at  the  same  time  as  it 
increases  wave  propagation  speed,  thereby  reducing  arrival 
time).  Therefore,  a stiffer  loading  behavior  was  introduced 
by  halving  the  parameter  W in  all  layers.  This  has  the  effect 
of  halving  the  volumetric  hysteresis,  or  locking  strain,  and 
increases  the  loading  moduli  of  the  material  while  leaving  the 
unloading  moduli  unchanged.  The  CIST  computation  was  repeated, 
thus  completing  the  first  cycle  of  the  iteration  procedure. 

Although  the  results  of  the  second  computation  were  much 
improved,  it  was  clear  that  further  improvement  could  be  obtained 
by  increasing  the  unloading  wave  speeds  in  the  model  (and, 
therefore,  the  parameters  K and  G).  This  was  done  and  several 
iteration  loops  were  completed  in  which  quantitative  variations 
of  K,  G and  W in  the  various  layers  were  made.  This  is,  of 
course,  a trial  and  error  process.  The  last  computation  of  the 
series  resulted  in  the  final  CIST  based  material  model.  The 
material  parameters  for  this  model  are  listed  in  Table  II  and 
the  behavior  of  the  model  is  shown  in  Figs,  (lb-25).  The 
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TABLE  1 I 


CAP  MODEL  BASED  ON  ITERATION  WITH  CIST 


MODEL 

PARAMETER 

DRY 

CLAY 

WET 

CLAY 

LIME  TONE 

SHALE 

COAL 

K (ks  i ) 

72  5 . 

1 020  . 

2 390. 

1160. 

580. 

G ( k s i ) 

2. 9 

6 . 09 

1440. 

247  . 

3 48  . 

A(ksi) 

.0334 

.0  16 

58.0 

55.2 

33.8 

- 1 

B ( k s i ) 

0. 

0 . 

. 00966 

. 00966 

. 00966 

' 

C ( k k i ) 

0. 

0. 

57  . 1 

5 5.1 

3 3.4 

- 1 

D ( k s 1 ) 

.483 

3.45 

. 69 

. 1 1 

. 586 

W 

1 

.054 

.0011 

. 0001  3 

. 0029 

| 

.ool  ! 

I R 

2.3 

2 . 5 

3.  5 

5 . 0 

3.8 

j 

Dens  lty(p.c. f. ) 

124  . 

128. 

167. 

I 

155. 

O' 

^3- 

oc 

AXIAL  STRESS  <KSI> 

0-00  *40  - 80  1-20  1*60  2-00  2*40 


UNIAXIAL  STRAIN  TEST  HARD  PAN  DRY  CLAY 
— * — cist-wkl 

LAI-  DATA 


FIG.  16 


COMPARISON  OF  LAB 


DATA  AND  CIST  BASED  MODEL 
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UNIAX  STRESS  PATH  HARD  PAN  DRY  CLAY 

— »■  - C 1ST -Mi DEL 
UW-  DATA 
FAILURE 


FIG.  17 


COMPARISON 


OF  LAB 


DATA  AND  CIST  BASED  MODEL 


AXIAL  STRESS  <KSI> 

0*00  -40  >80  1*20  1*60  2*00  2*40 


UNIAXIAL  STRAIN  TEST  HARD  PAN  WET  CLAY 

► CXST-MiDEL 

— ► — LAB  • BATA 


FIG.  18  COMPARISON  OF  LAB  DATA  AND  CIST  BASED  MODEL 
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FIG.  20  COMPARISON  OF  LAB  DATA  AND  CIST  BASED  MODEL 
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SIG1-SIG3  «SI) 

0-00  1-20  2-40  3-60  4-00  6-00  7*20 


UNIAX  STRESS  PATH  HARD  PAN  LIMESTONE 

► CIST -HA DEL 

LAI.  DATA 
FAILURE 


FIG.  21  COMPARISON  OF  LAB  DATA  AND  CIST  BASED  MODEL 
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0*00  -40  *80  1>20  1*60  2*00  2*40 


FIG.  22  COMPARISON  OF  LAB  DATA  AND  CIST  BASED  MODEL 
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FIG.  23 
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FIG. 25  COMPARISON  OF  LAB  DATA  AND  CIST  BASED  MODEL 
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computational  results  obtained  with  this  CIST  based  model  are 
shown  in  Figs.  (26-31).  They  are  labeled  "CALC.  WITH  C.  1 ST 
BASED  PROPS"  and  are  compared  with  the  corresponding  CIST 
measurements.  For  convenience  the  results  of  the  first  CIST 
computation  are  also  shown  in  these  figures  and  are  lableed 
"CALC.  WITH  LAB  BASED  PROPS". 

Figure  (26-31)  show  that  the  agreement  between  the  CIST 
measurements  and  the  final  CIST  computation  is  good.  The 
computed  magnitudes  and  waveforms  of  the  horizontal  velocity 
time -histories  throughout  the  region  of  the  CIST  event  compare 
well  with  the  overall  behavior  measured  in  CIST.  The  agreement 
of  the  CIST  measurements  to  the  computational  results  obtained 
with  the  LAB  based  model  is  of  course  not  as  good,  but  is  still 
substantial;  it  is  possible  that  the  LAB  based  model  for  the 
HARD  PAN  site  would  be  adequate  for  many  ground  shock  applications. 

It  is  important  to  note  that  the  CIST  based  model  behavior 
as  shown  in  Figs.  (16-25)  is  not  overwhelmingly  different  from 
the  original  laboratory  data.  It  is  reassuring  that,  at  least 
for  this  site,  a model  has  been  obtained  that  can  adequately 
represent  both  laboratory  and  in  situ  data. 

It  should  be  pointed  out  that,  although  all  of  the  CIST  results 
presented  here  involve  horizontal  motions,  measurements  and 
computations  of  vertical  ground  motions  are  also  available.  As 
one  might  expect,  however,  these  do  not  agree  well  at  all  because 
the  vertical  motions  in  CIST  are  extremely  sensitive  to  non- 
uniformity in  the  layering  and  to  material  non -homogeneity 
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FIG.  26 


HORIZONTAL-  COMPARISON  OF  MEASURED 
CALCULATED  CIST  WAVEFORMS  AT  Z = 45 
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FIG.  27  HORIZONTAL-  COMPARISON  OF  MEASURED  AND 
CALCULATED  CIST  WAVEFORMS  AT  Z = 35  FT 
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CALC.  WITH  CIST  BASEO  PROPS 
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FIG.  28  HORIZONTAL-  COMPARISON  OF  MEASURED  AND 
CALCULATED  CIST  WAVEFORMS  AT  Z = 35  FT 
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FIG.  29  HORIZONTAL-  COMPARISON  OF  MEASURED  AND 
CALCULATED  CIST  WAVEFORMS  AT  Z =29.5  FT 
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FIG.  30  HORIZONTAL-  COMPARISON  OF  MEASURED  AND 
CALCULATED  CIST  WAVEFORMS  AT  Z = 28 5 FT 
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FIG. 


* 


31  HORIZONTAL-  COMPARISON  OF  MEASURED  AND 
CALCULATED  CIST  WAVEFORMS  AT  Z = 18  FT 
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within  the  layers. 


Because  most  of  the  non -uniformities  are 


always  unknown  and,  even  if  known,  certainly  can  not  be 
included  in  two  dimensional axisymmetric  CIST  computations, 
it  seems  futile  to  hope  that  they  can  ever  be  represented, in  a 
CIST  computation  of  the  type  performed  in  this  study. 


48 


CONCLUSIONS  AND  RECOMMENDATIONS 


V 1 

By  means  of  an  iteration  procedure  in  which  a cap  model 
fitted  to  laboratory  data  was  altered  in  such  a way  as  to  lead 
to  improved  results  in  the  computation  of  CIST  12,  a model 
has  been  obtained  for  the  HARD  PAN  site  which  agrees  reasonably 
well  with  both  laboratory  and  in  situ  measurements.  This  model 
is  now  being  used  in  a number  of  ground  shock  and  structure- 
medium  interaction  computations  for  the  HARD  PAN  SITE. 

The  construction  of  the  final  cap  model  (which  is  consis- 
tent with  laboratory  as  well  as  CIST  data)  presented  little 
difficulty  in  this  study.  It  is  not  clear,  however,  whether 
this  is  a general  phenomenon  or  whether  it  is  merely  a fortunate 
consequence  of  the  behavior  of  the  HARD  PAN  site  in  particular. 
Therefore,  it  is  envisioned  that  a number  of  other  sites  at 
which  CIST  tests  have  been  performed  will  be  studied  using 
the  procedure  developed  in  this  investigation.  These  sites 
include  the  Nevada  Test  Sift"  (CIST  5,  6,  7)  the  MIDDLE  CUST 

wet  site  (CIST  9)  and  the  DICE  THROW  site  (CIST  15). 

The  present  study  has  shown  that  the  vertical  velocities 
observed  during  a CIST  test  are  too  sensitive  to  layer  irregu- 
larities to  be  useful  for  fitting  models.  For  this  reason, 
it  is  recommended  that  emphasis  should  be  placed  on  fielding 
additional  and/or  redundant  horizontal  velocity  gauges  in 
future  CIST  tests.  This  will  permit  one  to  obtain  estimates 
of  data  scatter  and  uncertainty,  and  to  establish  the  magnitude 
of  baseline  shift  which  can  be  expected  in  the  CIST  measurements. 
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Such  Information  would  enable  the  modeler  to  j u d <>  e better 
the  adequacy  of  any  particular  fit  and  to  evaluate  the 
accuracy  which  can  be  achieved  in  a ground  shock  computation 
using  the  model. 

Because  the  amount  of  useful  data  obtained  from  CIST 
is  limited  to  horizontal  motions  and  may  not  directly  relate 
to  deep  outrunning,  airblast  or  crater  induced  behavior,  the 
need  for  additional  information  is  clear.  Especially  when 
the  site  is  anisotropic,  tests  which  excite  primarily  vertical 
motions,  e.g.,  DISK  HEST  or  above  ground  bursts,  are  important 
so  that  laboratory  data  (which  is  obtained  mainly  from  verti- 
cally cored  samples)  may  be  related  to  the  in  situ  behavior 
in  that  direction.  Hopefully,  a comparison  of  laboratory  and 
in  situ  behavior  will  become  available  for  the  CIST  configuration 
when  TERRA  TEK,  Inc.,  completes  tests  in  which  CIST-like  stress 
paths  will  be  applied  to  laboratory  specimens  and  the  results 
compared  to  CIST  measurements  and  calculations. 

As  far  as  the  HARD  PAN  site  in  particular  is  concerned, 
it  is  recommended  that  the  material  model  obtained  in  this 
study  be  used  in  a calculation  of  the  DISK  HEST  which  was 
conducted  as  one  of  the  events  in  the  HARD  PAN  series.  A 
better  assessment  of  the  current  HARD  PAN  model  may  then  be 
made  and,  if  required,  an  improved  model  could  be  constructed 
to  take  into  account  this  additional  in  situ  information. 
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APPEND  I X 


Brief  Descript,  ion  of  t he  LAYER  Code 

The  LAYER  code  is  an  explicit  two  dimensional  axtsymmetric 
dynamic  finite  difference  code  designed  for  the  efficient 
solution  of  problems  involving  wave  propagation  in  a layered 
medium.  The  code  can  accept  a wide  range  of  material  model 
behavior  for  each  layer  of  a horizontally  bedded  geometry 
and  can  be  modified  to  accept  arbitrary  initial  conditions  and 
traction  and/or  velocity  boundary  conditions.  The  use  of 
reasonably  "quiet"  transmitting  boundaries  allows  half-space 
problems  to  be  efficiently  solved.  The  code  neglects  material 
t ran  sport  (convect ion)ef  fects  . 

A typical  grid  layout  for  the  LAYER  code  is  shown  in 
Fig.  A- 1 . Each  of  the  layers  consists  of  a homogeneous  single 
phase  transversely  isotropic  material.  Non- slip  contact  between 
layers  is  assumed,  although  slip  and  seperation  effects  can  be 
modeled  by  introducing  an  additional  thin  layer  of  a material 
which  is  weak  in  shear  and/or  tension. 

As  can  be  seen  from  Fig.  A - 1 , various  types  of  boundary 
conditions  may  he  applied.  These  Include  applied  surface 
pressures  as  well  as  rigid,  roller  and  "transmitting"  boundaries. 
In  addition,  the  annulus  shown  in  Fig.  A-l  may  be  replaced 
by  a solid  cylinder  by  setting  r = 0.  This  automatically 
introduces  the  proper  symmetry  and  continuity  boundary  condi- 
tions along  the  center  line.  Plane  strain  problems  may  be 
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PRESSURE 


FIG  A- 1 


TYPICAL  LAYER  CODE  GRID  LAYOUT 
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solved  by  us  inf  a large  value  of  r . 

Because  transport,  or  convective,  effects  are  ignored  in 

LAYER,  it  is  possible  to  formulate  the  governing  continuum 

equations  so  that  the  stress  and  velocity  are  the  only  unknown 

field  variables.  At  any  point  (r,  z ) , let  a , o , or.  and  t 

r z u 

represent  the  radial,  axial,  tangential  and  shear  stress  com- 
ponents in  the  cylindrical  (r,  z)  coordinate  system  as  functions 
of  r,  z and  time  t and  let  u and  v represent  the  radial  and 
axial  velocity  components  at  (r,  z,  t).  The  governing  equations 
consist,  of  the  equations  of  motion 
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in  which  p is  the  local  density  and  g is  the  acceleration  of 


gravity,  together  with  the  constitutive  equations,  which  may 
be  represented  as 
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The  finite  difference  analog  of  the  above  equations  is 
constructed  for  LAYER  by  using  a staggered  grid.  Fig.  (A-l), 
and  by  taking  central  differences  for  all  spatial  derivatives, 
forward  time  derivatives  in  the  equations  of  motion,  and 
backward  time  derivatives  in  the  constitutive  equations. 
Therefore,  at  any  time,  the  equations  of  motion  are  first 
employed  to  determine  the  new  velocity  field  and  this  field 
is  then  used  to  construct  the  strain  rate  components  3u/3r, 

3 v / 3 z , u/r  and  ( 9 u ^ 3 z + 9v/9r)/2.  These  are  then  used  in 
the  material  subroutine  to  obtain  the  new  stresses. 

The  I. AYER  computational  grid  is  orthogonal  with  a constant 
radial  mesh  size,  Ar,  anil  a variable  vertical  mesh  size,  Az, 
which  changes  arbitrarily  with  depth.  Because  long  and  narrow 
grid  elements  can  sometimes  introduce  significant  numerical 
errors  into  a computation,  the  aspect  ratio,  A r / A z , for  each 
mesh  point  is  generally  kept  between  1/5  and  5.  Of  course, 
the  usual  Co u rant,  Friedrichs,  Lewy  stability  condition  governs 
the  choice  of  time  step.  At. 

There  are  several  additional  features  of  interest  in 
I.AYER.  in  order  to  improve  efficiency  in  solving  certain  types 
of  extended  ground  shock  problems,  LAYER  has  a restart  capa- 
bility as  well  as  an  option  to  shift  the  computational  grid  in 
space  (dropping  grid  cells  behind  and  acquiring  undisturbed 
new  cells  in  front).  As  far  as  output  is  concerned,  quantities 
such  as  stress  and  velocity  components,  stress  invariants,  or 
other  desired  vari  allies,  are  written  at  each  time  step  and  at 
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specified  points  of  interest  on  a binary  tape.  Further,  peak 
stresses,  velocities,  and  displacements  are  stored  for  all 
depths  at  given  ranges,  and  peak  displacements  are  stored  for 
all  ranges  at  specified  depths.  Plots  of  the  time  histories 
and  cross  plots  of  pairs  of  quantities  can  easily  be  obtained. 

The  LAYER  code  has  been  implemented  on  a number  of  computers, 
but  it  has  been  used  principally  on  Control  Data  Corporation 
equipment,  specifically  on  CDC  6600  and  7600  machines.  A 
large  problem  may  have  about  10,000  space  points  and  involve 
about  2000  time  steps.  Such  a run  typically  requires  about 
three  hours  of  central  processor  time  on  a CDC  6600  if  the  cap 
model  is  used. 
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Chief  of  Naval  Re  at?  arch 
ATTN:  Tech.  Lib. 

ATTN:  Code  464,  Jacob  L.  Warner 
ATTN:  Code  464,  Thomas  P.  Quinn 
ATTN:  Nicholas  Perron© 

Officer-  in- Charge 

Civil  Engineering  Laboratory 
ATTN:  Stan  Takahashi 
ATTN:  Tech.  Lib. 

Commander 

Naval  Electronic  Systems  Command 
ATTN:  P ME  i 1 7- 2 1 A 

Commander 

Naval  Facilities  Engineering  Command 
ATTN:  Code  04 B 
ATTN:  Tech.  Lib. 

ATTN:  Code  03 A 

Superintendent  (Code  1424) 

Naval  Postgraduate  School 

ATTN:  Code  2124,  Tech.  Rpts.  Librarian 

Director 

Naval  Research  Laboratory 

ATTN:  Code  2027,  Tech.  Lib. 

Commander 

Naval  Sea  Systems  Command 
ATTN.  ORD-91313,  Lib. 

Commander 

Naval  Ship  Engineering  Center 
ATTN.  Tech.  Lib 

Commander 

Naval  Ship  Rsch  & Development  Ctr. 

ATTN.  Code  L42-3,  Library 

Commander 

Naval  Ship  Rsch  *.  Development  Ctr. 

Underwater  Explosive  Research  Division 
ATTN:  Tech.  Lib. 

Commander 

Naval  Surface  Weapons  Center 

ATTN:  Code  WX-21,  Tech.  Lib. 

ATTN  Code  W A- 801.  Navy  Nut*.  Prgms.  Off. 

Commander 

Naval  Surface  Weapons  Center 
ATTN:  Tech.  Lib. 

President 

Naval  War  College 

ATTN:  Tech  Lib 

Commanding  Officer 

Naval  Weapons  Evaluation  Facility 
ATTN:  Tech  Lib 

Director 

Strategic  Systems  Protect  Office 
ATTN  NSP-43,  Tech  Lib 


AF  Geophysics  Laboratory,  AFSC 

ATTN;  SUOL,  AFCRL,  Rsch.  Lib. 

AF  Institute  of  Technology,  Al’ 

ATTN:  Library,  AF1T,  Bldg.  640,  Area  B 

AF  Weapons  Laboratory,  AFSC 

ATTN:  DEV,  Jimmie  L.  Bratton 

ATTN:  DES,  M.  A.  Plamondon 

ATTN:  Robert  Port 

ATTN:  DYT 

ATTN:  Robert  Kenny 

ATTN:  SUL 

Headquarters 

Air  Force  Systems  Command 
ATTN:  DLCAW 
ATTN:  Tech.  Lib. 

Commander 

Foreign  Technology  Division,  AFSC 
ATTN:  TD-BTA,  Library 

HQ  USAF/IN 

ATTN:  IN  AT  A 

HQ  USAF/PR 
ATTN:  PRE 

HQ  USAF/RD 

ATTN:  RDQPN 

Commander 

Rome  Air  Development  Center,  AFSC 
ATTN:  EMTLD,  Doe.  Lib. 

SAM  SO/M  N 

ATTN:  MMI1 

Commander  in  Chief 

Strategic  Air  Command 

ATTN:  NRI-STINFO,  Library 

ENERGY  RESEARCH  & DEVELOPMENT  ADMINISTRATION 

University  of  California 

Lawrence  Livermore  Laboratory 
ATTN:  Tech.  Info.,  Dept.  L-3 

Los  Alamos  Scientific  Laboratory 

ATTN:  Doc.  Con.  tor  Reports  Lib. 

ATTN:  Doc.  Con.  for  G.  R.  Spillman 
ATTN:  Doc.  Con.  for  R.  J.  Br  id  well 

Sandia  Laboratories 

Livermore  Lalxiratory 

ATTN:  Doc.  Con.  for  Tech.  Lib. 

Sandia  Laboratories 

ATTN  Doc.  Con.  for  3141,  Sandia  Rpt.  Coll. 

t'.S.  Energy  Rsch.  & Dev.  Admin 

Albuquerque  Operations  Office 

ATTN:  Doc.  Con.  for  Tech.  Lib. 

U.S.  Energy  Rsch  l Dev.  Admin 

Div  ision  of  Headquarters  Services 

ATTN:  Doc.  Con.  for  Clans.  Tech.  Lib. 
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DEPARTMENT  OF  DEFENSE  CON' TRACTOR*  (Continued) 


U.S.  Energy  Rsch.  & Dev.  Admin. 

Nevada  Operations  Office 

ATTN:  Doc.  Con.  for  Tech.  Lib. 

Union  Carbide  Corporation 

Holifield  National  Laboratory 

ATTN:  Civil  Def.  Res.  Proj. 

ATTN:  Doc.  Con.  for  Tech.  Lib. 

OTHER  GOVERNMENT  AGENCIES 

Department  of  the  Interior 

Bureau  of  Mines 

ATTN:  Tech.  Lib. 

Department  of  the  Interior 

U.S.  ecological  Survey 
ATTN:  J.  H.  Healy 
ATTN:  Cecil  B.  Raleigh 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

Aerospace  Corporation 

ATTN:  Tech.  Info.  Services 

Agbabian  Associates 

ATTN:  M.  Agbabian 

Applied  Theory,  Inc. 

2 cy  ATTN:  John  C.  Trulio 

Avco  Research  & Systems  Group 

ATTN:  Research  Library,  A-H30,  Rm  7201 

Battelle  Memorial  Institute 
ATTN:  Tech.  Lib. 

The  BDM  Corporation 
ATTN:  Tech.  Lib. 

The  Boeing  Company 

ATTN:  R.  M.  Schmidt 
ATTN:  Aerospace  Library 

California  Research  ti  Technology,  Inc. 

ATTN:  Tech.  Lib. 

ATTN:  Sheldon  Shuster 
ATTN:  Ken  Kreyenhagen 

Cal  span  Corporation 
ATTN:  Tech.  Lib. 

Civil/Nuclear  Systems  Corp. 

ATTN:  Robert  Crawford 

University  of  Dayton 

Industrial  Security  Super  KL-505 
ATTN  . Hal  lock  F.  Swift 

University  of  Denver 

Colorado  Seminary 

ATTN:  Sec.  Officer  for  J.  Wi  sot  ski 
ATTN:  Sec.  Officer  for  Tech.  Lib. 

EGIrGt  Inc. 

Albuquerque  Division 
ATTN:  Tech.  Lib. 


General  American  Transportation  Corp. 

General  American  Research  Division 
ATTN:  G.  L.  Neidhardt 

General  Electric  Company 
TEMPO- (’enter  for  Advanced  Studies 
ATTN:  DA  SI  AC 

IIT  Research  Institute 
ATTN:  Tech.  Lib. 

Institute  for  Defense  Analyses 

ATTN:  [DA,  Librarian,  Ruth  S.  Smith 

Kaman  AviDyne 

Division  of  Kaman  Sciences  Corp. 

ATTN:  Tech.  Lib. 

ATTN:  E.  S.  Criscione 

Kaman  Sciences  Corporation 
ATTN:  Library 

Lockheed  Missiles  & Space  Co.,  Inc. 

ATTN:  Tech.  Lib. 

Lockheed  Missiles  & Space  Company 
ATTN:  Tech.  Info.  Ctr. , D/Coll 
ATTN:  Tom  Geers,  Dept.  52-33,  Bldg.  205 

McDonnell  Douglas  Corporation 
ATTN:  Robert  W.  Halprin 

Merritt  Cases,  Incorporated 
ATTN:  J.  L.  Merritt 
ATTN:  Tech.  Lib. 

The  M it  re  Corporation 
ATTN:  Library 

Nathan  M.  Newmark 
Consulting  Engineering  Services 
ATTN:  Nathan  M.  Newmark 

Physics  International  Company 

ATTN:  Doc.  Con.  for  Tech  Lib. 

ATTN:  Doc.  Con.  for  Larry  A.  Behrmann 
ATTN  DM  ON  t<»r  Ini  \1  | 

ATTN:  Doc.  Con.  for  E.  T Moore 

ATTN:  Dm.  On  Imp  MmI  Wi 

ATTN:  Doc . Con  for  Charles  Cktdfrey 
ATTN:  Doc.  Con  for  Dennis  Orphal 

R & I)  Associates 

ATTN:  Tech.  Lib 
ATTN:  Harold  L.  Brink- 
ATTN:  Cyrus  P.  Knowles 
ATTN:  Bruce  llartenbaum 
ATTN:  Albert  L.  Latter 
ATTN:  Henry  Cooper 
ATTN:  J.  G.  Lewis 
ATTN:  William  B.  Wright,  Jr. 

ATTN:  Jerry  Carpenter 

Science  Appl  ications,  In 

ATTN:  D.  E.  Maxwell 
ATTN:  David  Bernstein 

Science  Applications,  Inc 
ATTN  Michael  McKay 
ATTN:  Tech.  Lib. 
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DEPARTMENT  OE  DK1KNSK  CONTRACTORS  (Continued) 


s icnt'i  Applications,  Inc. 
ATTN  K.  A.  Shunk 


Universal  Analytics,  Inc. 
ATTN:  E.  I.  Field 


Soothwr-i  Research  Institute 
ATTN:  Wilfred  E.  Baker 
ATTN:  A H.  Wenzel 

Stanford  Research  Institute 

ATTN;  George  R.  Abraham  son 
ATTN:  Carl  Pel*  r son 
ATTN:  Burt  R.  Gasten 

Systems,  Science  K Software,  Inc. 
ATTN:  Tech.  Lib 
ATTN  Ted  Cherry 
ATTN.  Thomas  L).  Riney 
ATTN:  Donald  R.  Grine 

Terra  Tek,  Inc-. 

ATTN:  Tech.  Lib. 

ATTN:  Sidney  Green 


URS  Research  Company 
ATTN:  Tech.  Lib. 

The  Eric  H.  Wang  Civil  Engineering  Rseh.  Fac. 
ATTN:  Neal  Baum 
ATTN:  Larry  Bickle 

Washington  State  University 
Administrative  Office 

ATTN:  Arthur  Miles  Hohorf  for  George  Duval 


Weidlinger  Assoc.  Consulting  Engineers 


ATTN 

ATTN 

ATTN 


Melvin  L.  Baron 
J.  W.  Wright 
Dr.  I.  S.  Sandler 


Weidlinger  Assoc.  Consulting  Engineers 
ATTN:  J.  Isenberg 


Tetra  Tech,  Inc. 

ATTN:  Tech  Lib. 

ATTN:  Li- San  Hwang 

TRW  Systems  Group 

ATTN:  D.  H.  Baer,  Rl-2136 
ATTN:  I.  E.  Alter,  Rl-1008 
ATTN:  lech.  Info.  Ctr.,  S-1930 

ATTN:  R.  K.  Plehuch,  Rl-2078 
2 cy  ATTN.  Peter  K.  Dai,  R 1-2170 

TRW  Sy  stems  Group 
San  Bernardino  Operations 

ATTN:  E.  Y.  Wong.  527/712 


Westinghouse  Electric  Company 
Marine  Division 

ATTN:  W.  A.  Volz 
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